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A STUDY OF THE POSITIONS AND VELOCITIES OF A SPACE STATION AND A 
FERRY VEHICLE DURING RENDEZVOUS AND RETURN 

By John M. Eggleston and Harold D. Beck 


SUMMARY 

A study is made of the families of nonthrusting 
ascent trajectories of a ferry vehicle during rendez- 
vous with an orbiting body, referred to as a space 
station . It is shown that these trajectories may also 
be interpreted as descent trajectories of the ferry from 
the station to the earth. The rendezvous trajectories 
start at the end of the boost period {assumed to 
be 60 miles) and terminate at the station. The 
equations of motion are derived and results are shown 
for two typical orbits of the station: a 300-mile cir- 
cular orbit and a 100- to 500-mile elliptical orbit. 
Trajectories are described in terms of a rotating co- 
ordinate system fixed in the station and launch con- 
ditions are tabulated in terms of nonrotating 
inertial coordinates. 

Boundaries are given in terms of launch {at time 
of booster burnout) and rendezvous conditions for 
the example cases. The considerations used to cal- 
culate these boundaries and the significance of some 
of the trends are discussed. 

INTRODUCTION 

In its essence t lie rendezvous problem may be 
formulated in terms of two vehicles, the space 
station and the ferry vehicle. The space station 
is in an established orbit. The ferry vehicle is to 
be launched from some point on the earth and 
placed in a trajectory which intersects the orbit of 
the space station at the same instant that the space 
station reaches this intersection point. An addi- 
tional requirement, that the relative velocities of 
the two vehicles be zero at interception, may or 
may not be imposed. If such a requirement were 
imposed, the zero relative velocity would ordi- 
narily be obtained just prior to interception by 
terminal thrust control. Although formulated in 
terms of the approach of the ferry to the station, 


it is evident that the related problems such as that 
of ejection of the ferry or some other mass from the 
space station so that it arrives above a certain 
point on the earth and the interception of orbital 
bodies do not involve any essential difference with 
rendezvous as described herein. 

One of the primary concerns associated with the 
rendezvous mission is that of achieving once-a-day 
rendezvous capability and still not spending a 
large percentage of the payload on fuel for ex- 
pensive space maneuvers. (See ref. 1.) If the 
thrusting maneuvers are to be kept small, the sta- 
tion must be in a favorable position in its orbit at 
the time the ferry launch site passes through or is 
in the near vicinity of the plane of the orbit. The 
problem may be resolved into two parts: (1) 
placing the station into a desirable orbit so that a 
favorable position of the station near the ferry 
launch site occurs once a day and (2) extending 
the time during which the station is in a favorable 
position to allow for delays in ferry launch time. 

The first of these problems has been studied in 
references 1 to 4. References 1 and 2 point out 
that, if the maximum latitude of the station is 
only slightly greater than the latitude of the ferry 
launch site, once-a-day out-of-plane launches may 
be made with only a slight reduction in mass ratio. 
References 3 and 4 point out that rendezvous- 
compatible orbits for the station may be obtained 
if the station period is integrally related to the 
rotational period of the earth. 

One approach to the second part of the problem 
is to avoid restriction of the ascent trajectory of 
the ferry to a minimum-energy or Hohmann ellipse 
but to investigate the variations in the ferry trajec- 
tory and in the relative positions of the ferry and 
station at the time of launch because of variations- 
in the launch conditions (velocity, flight-path 

1 


2 


TECHNICAL REPORT R— 87* — NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


angle, and time). The restrictions that rendez- TJ 
vous must still be achieved and achieved at a AV 
reasonable cost in mass ratio are then boundary 
conditions of the problem. This paper is con- V 
corned with the variations in the trajectories and 
launch conditions which will lead to a successful V r 
rendezvous. V t 

When the various rendezvous trajectories were Ve 

considered, no attempt was made to eliminate 
those launch conditions which do not lend them- V r 

selves to “safe-launch” trajectories or to the limi- 
tations of specific boosters, although such consid- A W 
erations should be taken into account for any 
specific rendezvous mission. W 


X,Y,Z 

SYMBOLS 


Any consistent set of units may be used for the 


analysis 

presented herein. In the examples given, 


it is assumedthat 


<7 

e =32.17 feet per second per second 
r e =3,960 statute miles 

z, y,z 


1 statute mile = 5,280 feet 
1 foot = 0.3048 meter 

a 

A 

seinimajor axis of elliptic orbit of 



space station 


C 

constant angular momentum per unit 
mass 

E 

total energy 


F t 

generalized external forces (lift, drag, 
thrust) 

y 

g 

acceleration due to gravity 

e 

9e 

gravitational constant at earth's sur- 
face 

r 

h 

altitude above earth's surface 


k 

unit vectors in direction of X-, F-, and 


Z-axes, respectively 

Vf 

I sp 

specific impulse 


L 

Lagrangian, defined as T~U 

3 

m 

mass of ferry vehicle 


<Li 

generalized coordinates in Lagrangian 
operation, i= 1, 2, 3 


R 

distance from space station to ferry 



vehicle 

e. 

r 

radial distance measured from earth's 
center of gravity 


T Z , Ty, 

T z components of vehicle thrust acting in 
x-, y-, and ^-directions, respectively 

A0 

T 

kinetic energy 


t 

time 



potential energy 

relative velocity of ferry with respect 
to station 

total velocity of ferry vehicle as ob- 
served from center of earth 
radial component of V 
tangential component of V 
component of V t parallel to plane of 
station orbit- 

component of V t perpendicular to 
plane of station orbit 
change in the ferry vehicle weight due 
to fuel expenditure 
weight of ferry vehicle 
axes system with center fixed in space 
station. The Z-axis is perpendicu- 
lar to plane of orbit, the X- and F~ 
axes lie in plane of orbit with F 
always pointing away from center 
of earth. 

three orthogonal components of rec- 
tangular coordinate system 
angle measured from X-axis to projec- 
tion of relative velocity AF on the 

x f y plane, tan -1 ^ 

angle between relative velocity AV 

and the x,y plane, tan" 1 — ~ — _ 

yx 2 +?/ 2 

flight-path angle between local horizon 
and V , positive upward 
eccentricity of space station orbit 
separation angle between ferry and 
orbital plane of station as measured 
from center of earth, positive when 
station is on right side of ferry 
angle between orbital plane of space 
station and velocity component V t 
angular position of space station as 
measured in plane of its orbit about 
center of earth. (In an elliptic orbit, 
0=0 at perigee; in a circular orbit, 
0=0 at rendezvous.) 
angular position of ferry, measured in 
plane of station's orbit from posi- 
tion 0=0 

separation angle between ferry and 
station as measured in plane of sta- 
tion orbit from center of earth, posi- 
tive when station is ahead of ferry 


A STUDY OF THE POSITIONS AND VELOCITIES OF A SPACE STATION AND A FERRY VEHICLE 3 


ft 

co 

Subscripts: 

a 

c 

e 

1 

o 

P 

s 

L 

max 


a specific value of time 
total angular velocity of space station 
constant chosen to be exactly equal to 
0 when space station is in exactly a 
circular orbit 

apogee position 

circular 

earth 

ferry vehicle 

initial values at a specific position 
perigee position 
space station 

launch values where launch occurs at 
end of boost period 
maximum 


Derivatives with respect to time are denoted 

do 

with dots over the variables. For example, ^-=9. 
Vectors are denoted by bold face letters. 


ANALYSIS 

MISSION CONCEPT 

The rendezvous mission is herein separated into 
two major phases: the rendezvous with another 
vehicle in an established orbit, and the return 
from that orbit back to earth. If the problem is 
restricted to a coasting phase outside of the sensi- 
ble atmosphere of the earth, both problems can 
be studied simultaneously by investigating the 
trajectories in both positive and negative time. 
This approach will be the one used in this analysis. 

The study is based on the assumption that the 
space station is in an arbitrary orbit and that the 
boost vehicle for the ferry is sufficiently powerful 
to launch the ferry into a variety of coasting tra- 
jectories, one of which is exactly right for the 
rendezvous based on the position of the station 
at the time of launch (booster burnout). The 
family of trajectories to be considered is illustrated 
in figure 1. As one deviates from a Hohmann 
ellipse (shown on the left-hand side of fig. 1), 
there is a family of transfer trajectories with 
apogee greater than the orbital altitude of the 
station or with perigee within the sensible atmos- 
phere. On such trajectories rendezvous is to be 
achieved either during the ascending or during 
the descending phase of the transfer trajectory. 

The problem investigated is essentially that of 


determining the range of conditions over which a 
ferry can be launched and can still follow a tra- 
jectory which terminates at the space station. 
Rather than launch the ferry over a range of con- 
ditions and eliminate those trajectories that do 
not lead to rendezvous, the problem is treated in 
reverse. All calculations start with the ferry 
adjacent to the station at time zero. The ferry is 
given an incremental velocity component relative 
to the station and the equations of motion are 
solved in negative time for the rendezvous phase 
and in positive time for the return phase. (It 
will be shown in a subsequent section that with 
the proper interpretation, both results can be 
obtained by considering only one problem — either 
rendezvous or return.) In all cases the velocity 
and position of the ferry are calculated until the 
ferry vehicle reaches some reference altitude above 
the earth, taken to be 60 statute miles or until a 
minimum altitude is reached. The conditions at 
the 60-mile altitude are considered launch condi- 
tions for the rendezvous mission or entry conditions 
for the return mission. 

The term launch is used to refer to the end of 
the boost phase and the start of the coasting phase 
on the ascent trajectory. Although current boost 
vehicles cannot, in general, achieve near-orbital 
velocities at 60 miles, future boost vehicles ap- 
proach this condition. Reference 5 describes the 
desirability of achieving booster burnout at an 
altitude (above the atmosphere) as low as possible 
and coasting up to orbital altitude where a small 
additional velocity is applied for injection into 
orbit. 

In this study it is not required that the relative 
velocity between the two vehicles be zero at the 
time of rendezvous since it is assumed that thrust 
control will be employed during the terminal phase 
to insure that the relative velocity and displace- 
ment are brought to zero simultaneously. How- 
ever, in establishing the ground rules for this 
study, it is considered that this velocity correction 
should be of reasonable magnitude and should 
require a discharge of mass that would not exceed 
10 or 15 percent of the mass of the launched ferry 
vehicle. 1 

It is further assumed that the trajectories of the 
two bodies (the station and ferry vehicle) do not 

1 A transfer ellipse having a perigee at 60 miles and an apogee at 300 miles 
above the earth will require a 366 feet per second increase in velocity for 
injection into a 300-mile circular orbit. With I» v of 250 seconds, this value 
represents a minimum reduction in mass of about 4.0 percent. 
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Figure 1. — Three typical trajectories of a ferry during a rendezvous with an orbiting space station. Upper figures 
show the trajectories as seen by an inertial observer; lower figures show the trajectories as seen by an observer in 
the space station. Station coordinates are X and Y; AV is the relative velocity of the ferry with respect to the 
station at the time of rendezvous. 


have to be in the same plane at the time of launch 
but their planes should intersect prior to or at the 
point of interception. The alternative (where the 
ferry is injected into an intersecting orbit) would 
make the injection guidance more difficult and the 
out-of-plane correction after injection more expen- 
sive. Therefore, in this study the relative velocity 
vector of the ferry vehicle at interception does not 
necessarily lie in the plane of the station's orbit. 
The implication in these last two assumptions is 
that terminal and probably inidcourse guidance 
will be available to produce minor modifications 
to the nominal coasting trajectories obtained 
herein. 

EQUATIONS OF MOTION 

The equations of motion employed describe the 
motions of a mass, the ferry, moving in a radially 
symmetric gravitational field whose center is the 
center of the earth. The motions of the mass are 


measured in a rotating set of coordinates whose 
origin is located at the center of the space station. 
The space station moves about the center of the 
gravitational field in a Kepleriari trajectory at an 
angular velocity of 0 which is a constant w only 
if the trajectory is a circle. The coordinate 
system moving with the station is composed of a 
rectilinear set of axes with the Z-axis normal to 
the orbital plane. The X - and F-axes are rotated 
about the Z-axis at the same angular velocity 
0 so that the F-axis always points away from the 
center of the earth. A schematic drawing of the 
coordinates of the problem are showm for the 
station in a circular orbit (fig. 2(a)) and in an 
elliptic orbit (fig. 2(b)). 

The equations of motion of a mass with respect 
to a coordinate system fixed in a space station 
moving in a noncircular orbit (the more general 
case) are derived in appendix A. The variations 
of the station coordinates r, and 0 are then 
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(a) Station in circular orbit. 

(b) Station in elliptic orbit. 

Figure 2. — Coordinates employed in describing the mo- 
tions of mass m in the space-station coordinate system. 

described by the Keplerian equations of motion 
for which the solutions are well-known. For the 
special case where the space station is in a circular 
orbit and r s and 0 are constants, the equations 
of motion reduce to a form given in appendix 
B. Also given in appendix B are some approxi- 
mate solutions to these equations which may be 
used when the distance from the mass to the 
origin of the coordinate system is small compared 
with r f . These approximate solutions were re- 
ported in references 6 and 7. Since neither 
reference 6 nor reference 7 gives the derivation, it 


has been included in appendixes A and B. 

The inertial position and velocity of the ferry is 
also important for various purposes such as the 
calculation of launch conditions. Transforma- 
tions for the position and velocity data from the 
rotating orbiting frame of axes to the nonrotating 
earth-fixed axes are given in appendix C. 

DESCRIPTION OF MOTION 

Since the motion of a body with respect to a 
rotating coordinate system moving in a central 
gravity field has certain unique characteristics, 
some description of that motion is given. Con- 
sider at time zero a mass located at the center 
of the rotating coordinate system where the 
system is moving in a circular orbit. If that 
mass has an initial relative velocity AV in the 
plane of the orbit, the mass will at subsequent 
times move out from the origin on a curved path. 
This condition is illustrated in figure 3 for the 
case of the coordinate system moving counter- 
clockwise in a circular orbit of radius r,=4,260 



Figure 3. — Trajectories of a number of point masses 
ejected from the center of the rotating coordinate system 
at J=0 (0 = 0), each with a total relative velocity of 
10 feet per second, but with different velocity compo- 
nents. The solid lines are discrete trajectories; the 
dashed lines are the contours of the positions of the 
masses at subsequent positions of the coordinate system. 
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miles out from the center of the earth (h s = 300 
miles). In figure 3 trajectories are shown for 
the mass with several different combinations of 
initial velocity components. In each case the 
total relative velocity is 10 feet per second. The 
angular velocity of the coordinate system to is 
constant so that the subsequent position of the 
mass is defined at any time t when the coordinate 
system has moved c ct=Q° around the earth. 
Dotted lines have been used to join the positions 
at certain values of 0 up to 0 — 90°. It can be 
seen that, if a mass has a relative velocity radially 
outward from the earth (along the F-axis so that 
x 0 =Q, 2 /o= 10 ), then, as its radial distance from 
the center of the earth becomes larger, the angular 
velocity of the mass becomes less than that of 
the coordinate system, and the mass moves up 
and behind the coordinate system. Likewise a 
mass having a relative velocity along the positive 
X-axis (x o —\0, y 0 = 0) will initially have a lower 
angular velocity than the coordinate system and 
will move back and below the coordinate system. 
As the mass moves closer to the earth, the angular 
velocity increases and the mass moves ahead of 
the coordinate system. For other combinations 
of x and y the mass will follow similar trajectories. 
The locus of points on these trajectories at equal 
values of time (or 0), shown by the dashed lines 
in figure 3, is first a circular pattern and later an 
elliptical pattern about the coordinate center. 

It would appear that the trajectories of figure 
3 are antisymmetric. For small relative velocities 
and values of 0 = c ot, this is approximately true. 
However, the spherically symmetric gravity field 
changes this condition; thus, the motions are 
actually antisymmetric about a curvilinear coor- 
dinate system defined by the F-axis and the line 
of constant r s Q as shown in figure 4. 

The effect of this spherical gravity field on the 
motions of a mass is shown in figure 5 for two 
special cases. In figure 5(a) the mass has an 
initial velocity of 400 feet per second directed 
along the F-axis of the coordinate system. In 
figure 5(b) the mass has the same velocity but the 
velocity is directed along the positive A"-axis of 
the coordinate system. Shown on the left-hand 
side of figure 5 are the relative motions of the mass 
and the coordinate system as seen by an inertial 
observer. On the right-hand side the same tra- 
jectory is shown as viewed by an observer in the 
center of the coordinate system. In order to 



Ei CURE 4. — Rectilinear and curvilinear coordinates. 1 


simplify the definition of the direction that AV 
makes with the coordinate system, the angles a and 
ft are introduced and are defined as 

a=t-an — 
x 


Tims, /3 is the angle that AV makes with the 
x,y plane of the coordinate system and a. is the 
angle that the in-plane component of AV makes 
with the /-axis. If AV lies entirely in the plane 
of the orbit of the coordinate system (z=z— 0), 
then a is simply the angle that Al 7 makes with 
the j’-axis. 

In figure 5(a), the value of a is 90° and the 
trajectory starts at 0=0°. The trajectory of the 
mass is an ellipse with apogee at 0 = 90° and 
perigee at 0 — 270°. When 0 is 360°, the mass 
has returned again to the center of the coordinate 
system. In this special case, the mass had the 
same orbital period as the coordinate system. 
This condition is the result of assuming a spherical 
gravity field. If the oblateness of the earth were 
taken into account, the orbital periods would not 
be the same. In a spherical gravity field this con- 
dition is always true for any relative velocity of 
the mass, when the angle a is equal to 90° or 270°. 

Figure 5(b) shows the case where a=0°. The 
trajectory starts with 0=0° and the inertial 
velocity of the mass 400 feet per second less than 
that of the coordinate system. Because of this 
smaller inertial velocity, the centrifugal force of 
the mass is less than that of the coordinate system 
and the mass initially falls behind and below. 
The period of the mass is shorter than that of the 
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9 = 90® 





(b) x , miles 


(a) AF—400 feet per second; a = 90°. 

(b) AT r =400 feet per second : « = 0°. 

Figure 5. — Two typical trajectories of a mass which at £=0 (0=0) was at the center of the coordinate system and had 
a relative velocity AV, which made an angle a with the X-axis. Trajectories are shown as viewed by an inertial 
observer (on the left) and an observer in the center of the coordinate system (on the right). 


coordinate system in this case; therefore, after the 
coordinate system has made one complete revo- 
lution about the earth, the mass is roughly 1,200 
miles ahead of the F-axis and 180 miles below the 
X-axis. Because of the curvature of the earth, 
however, both the mass and the coordinate system 
are at approximately the same altitude of 300 
miles. The curvilinear coordinate r s Q is shown 
in figure 5(b) as the line of constant h s which is, 
in this case, 300 miles above the mean earth. 

COMPARISON OF MOTIONS OBTAINED FROM EXACT AND 
APPROXIMATE EQUATIONS 

Since the equations of motion of a mass in terms 
of the rotating coordinate system are amenable to 
approximate solutions, it may be desirable at this 
point to examine the motions as computed from 
the approximate equations and to consider whether 

597S70 — -61 2 


motions so calculated are satisfactory for use in 
the problem under consideration. 

In appendix B, approximate solutions are ob- 
tained for the case where the coordinate system 
is in a circular orbit. The approximate equations 
were obtained by expanding the gravity term into 
a Taylor series (in powers of x> y , and z) and drop- 
ping all terms higher than the first order. The 
effect of this approximation on the calculated 
motions of the mass is illustrated in figure 6. At 
tf=0 the mass was at the center of the coordinate 
system and had a relative velocity of 200 feet per 
second in the direction «=45°. Subsequent tra- 
jectories were computed in both positive and 
negative time by using both the approximate and 
the exact equations of appendix B. In the upper 
half of figure 6 is shown the variation of y with x 
of the trajectory in terms of the rotating coordi- 
nate system. Shown in the lower half of figure 6 
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miles 



© = w/ , deg 


Figure 6. — Comparison of two trajectories obtained by solving exact and approximate equations of motion in both 
positive and negative time. At 2 — 0 both the approximate and the exact equations had the same initial conditions: 
AF=200 feet per second, a=45°, 0 = 0°. 


s the variation of the altitude of the mass desig- 
nated as hf , the distance above the mean surface 
of the earth. It may be seen that, for values of 
|e|^90°, the differences between the approximate 
and exact solutions are fairly small. By the time 
the coordinate system has made one complete 
revolution about the earth, it may be seen that 
the differences between the two solutions become 
very large. Within the first order of magnitude, 
this error is directly a function of 0 and not of x 
or y . The trajectory of the mass as computed by 
the approximate equations appears to be the result 
of a gravity field parallel to the F-axis. The 
approximation that the gravity field is parallel 
while the earth and actual gravitational field are 


curved results in an apparent increase in h f with 
each orbital revolution. With these approximate 
equations, errors in velocity and acceleration are 
obtained similar to those obtained in the position 
of the mass. 

From this cursory examination, it would appear 
that motions obtained from the approximate 
equations are usually adequate when the mass is 
in close proximity to the coordinate system and 
0 is less than, for example, 90°. However, for 
the case under consideration 0 and the separation 
distance may become relatively large, and it would 
appear that motions computed from the approxi- 
mate equations would not be sufficiently accurate. 
Therefore, only the exact equations are used in 
the subsequent sections of this paper. 
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SOLUTIONS OBTAINED IN POSITIVE AND NEGATIVE TIME 

In order to obtain a rendezvous between a 
space station already in orbit and a ferry vehicle 
which is launched from the earth, the ferry vehicle 
must be put onto a trajectory which intersects or 
nearly intersects the orbital track of the space 
station. Likewise, when the ferry leaves the space 
station for the return trip to some desired landing 
site on earth, it must be put onto a trajectory 
which intersects the earth’s atmosphere at some 
predetermined position. Both trajectories may 
be computed by placing the mass, which in this 
case is the ferry, at the center of the station’s 
coordinate system and running the problem in 
negative or positive time. Figure 7 illustrates 
two such trajectories calculated in both positive 
and negative time. 

One trajectory, shown in figure 7 with a solid 
line, represents the trajectory of the ferry vehicle 


for the case where AV was 200 feet per second and 
a was —135° at 2=0. It may be seen that for both 
positive and negative time the trajectory of the 
ferry was at a greater orbital altitude than that of 
the space station. Ticks have been placed on the 
trajectory for every integral multiple of 9 of 90°, 
that is, every quarter revolution of the station 
about the earth. Although not shown, a tra- 
jectory obtained for a=135° also showed the 
same characteristic of moving away from the 
earth in either positive or negative time. 

The trajectory shown with the dashed lines in 
figure 7 is a case where again AV=200 feet per 
second but a— 45° at 2=0. In this case the 
ferry reached a perigee altitude of 200 miles 
above the surface of the earth. This same 
perigee altitude was obtained for trajectories 
calculated in both negative and positive time and 
for a= —45° as well as for a= 45°. Although' this 



a = -135° 
a = 45° 


( t 



Figure 7. — Variation of two ferry vehicle trajectories computed in positive and negative time. Trajectories are shown 
in terms of the rotating coordinate system {x,y) of a space station in a circular orbit and in terms of altitude measured 
from the surface of the earth. At t~0 the ferry is at the center of the station coordinate system and has a relative 
velocity of 200 feet per second which makes an angle a with the Y-axis. Ticks are shown at every integral multiple 
of 0 = 90°. 
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altitude is still above that usually used to define 
the upper limit of the earth's atmosphere, the 
important point is that, when |a| was 45°, the 
ferry moved toward the earth, and, when \a\ was 
135°, the ferry moved away from the earth. 

A general rule may be stated here for the case 
of a station in a circular orbit: For either ren- 
dezvous or return trajectories and for a given 
relative velocity AV at the time of intersection of 
the station and ferry, the trajectory leading to 
the smallest perigee will have a relative velocity 
AV which makes an angle a=0° with respect to 
the X-axis (that is, the case where the trajectory 
of the ferry is a Hohmann ellipse). As \a\ in- 
creases, the perigee altitude will increase until at 
|a| = 180° the perigee altitude of the ferry is equal 
to the orbital altitude of the station. Thus, if the 
ferry is to go from the station to the proximity of 
the earth, or vice versa, with a near-minimum 
amount of energy, the relative velocity at the time 
of intersection will be in the first or fourth quad- 
rants (— 90°<a<90°). If AV is large enough 
and lies in the second or third quadrants, close 
proximity to the earth can be achieved but the 
inertial velocity will always be larger than if AV 
were in the first or fourth quadrants. 

In the case of a= — 135° in figure 7, it may be 
seen that the tnijectorj' had a perigee only a few 
miles closer to the earth than the altitude of the 
station. In the case of a=45°, perigee was 100 
miles closer to the earth. The same would be true 
for 135° and a=— 45°. Therefore, most of 
the trajectories studied in the remainder of this 
paper will be those where a lies in the first or 
fourth quadrants (|a|^90°). 

RESULTS 

By using the equations derived in appendixes 
A, B, and C, numerical results have been obtained 
for two particular orbits of a space station. In 
the first case, the station was assumed to be in 
a circular orbit 300 miles above the earth. In the 
second case, the station was assumed to be in an 
elliptic orbit having a perigee of 100 miles and 
an apogee of 500 miles so that the nominal orbital 
altitude was 300 miles. Only results for the 
rendezvous mission are shown and discussed. 
However, the associated case of the return trajec- 
tory of the ferry from the station to the earth 
(obtained by solving the equations of motion in 


positive time) are symmetric reflections (about the 
F-axis) of the rendezvous trajectories. Thus the 
rendezvous trajectories to be shown can also be 
interpreted as return trajectories by simply revers- 
ing the direction of AV and the direction of the 
station rotation. This interpretation is illustrated 
in figure 8. With this interpretation the launch 
conditions discussed for the rendezvous mission 

Direction of 
station rotation 



Direction of 
station rotation 



Figure 8. — Rendezvous aud return trajectories. 
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can be also treated as the entry conditions for the 
return mission. 

STATION IN A CIRCULAR ORBIT 

The trajectories leading to rendezvous of a ferry 
with the space station in a 300-mile circular orbit 
were obtained by placing the ferry at the center 
of the station’s coordinate system at £=0 and 
solving the equations of motion in negative time. 
Five cases of relative velocity &V were studied: 
200, 400, 600, 800, and 1,000 feet per second. The 
angle a was varied from —90° to 90°. The results 
of these trajectories are shown in figures 9 to 13. 

Coplanar trajectories. — The first trajectories 
presented are those obtained for the case where 
the entire motion of the ferry vehicle takes place 
in the plane of the station’s orbit. Shown in figure 
9 are a number of typical trajectories where the 
motion is described in terms of the x,y coordinate 
system of the station. It was assumed in these 
trajectory calculations that booster burnout oc- 
curred just above the sensible atmosphere of the 


earth, an altitude which was chosen to be 60 miles. 
Those trajectories which intersect the 60-mile alti- 
tude were considered to be possible rendezvous 
trajectories. Ticks are placed on the trajectories 
at this point. The results for the case of AF=200 
feet per second are not shown since none of the 
trajectories constitute possible rendezvous trajec- 
tories from a booster burnout altitude of 60 miles. 
From orbital mechanics calculations it may be 
shown that the minimum injection velocity (that 
is, closing velocity) is approximately 360 feet per 
second for the in-plane case considered. In figure 
9(a), AF=400 feet per second and trajectories 
having an |a|^20° are shown to be possible ren- 
dezvous trajectories. For |a|^30°, rendezvous 
trajectories from an altitude of 60 miles were not 
possible. Figure 9(a) illustrates that at booster 
burnout the ferry can be anywhere from 335 to 
415 miles behind the space station, a spread of 80 
miles. Since the station is traveling at a rate of 
approximately 5 miles per second over the surface 
of the earth, this distance represents a spread in 



(a) aF= 400 feet per second. 

Figure 9. — Variation of y with x for a series of trajectories of a ferry vehicle during rendezvous with a space station in a 
300-mile circular orbit. The ferry intersects the space station with a velocity AV and an angle of closure a. All 
trajectories art* in the orbital plane of the station. Ticks are used to denote those trajectories which start at a 60-mile 
altitude. 
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(b) A T r = 600 feet per second. 
Figure 9. — Continued. 


launch time of 16 seconds. Thus, as the angle 
of closure a opens up, so does the range of launch 
conditions and launch times. A similar analysis 
may be made from figures 9(b), 9(c), and 9(d) 
where AT 7 equals 600, 800, and 1,000 feet per 
second, respectively. 

Also shown in these figures is a tabulation of 
the launch conditions for each of the possible 
rendezvous trajectories. For each of these trajec- 
tories is listed the inertial velocity and flight-path 
angle of the ferry at launch, and the positional 
angle of the station at the time of launch 0^. By 
definition, 0 at rendezvous is equal to zero. Again 
referring to figure 9(a), it may be seen from the 
tabulated data that from the time of booster 


burnout to the time of rendezvous the space station 
travels between 147° and 162° around the earth. 
These values of launch conditions were obtained 
by simultaneously computing the inertial velocities 
by the method described in appendix C and by 
interpolating between the numerical data points 
to obtain the values at approximately 60 miles. 

The trajectories were initially obtained at in- 
crements of a of 30°. In order to investigate the 
trajectories that lay between possible and not 
possible rendezvous trajectories, additional inter- 
mediate cases were also calculated in 10° incre- 
ments. Thus, the pattern of values of a investi- 
gated for any one case may be irregular. 



miles 
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In figure 10 the variation of range rate with 
range is shown for the same trajectories that were 
treated in figure 9. Ticks are again put on the end 
of each trajectory at the point where it intersected 
an altitude of 60 miles. It is of interest to note 
the wide variation of range rate with range between 
the trajectories which led to rendezvous. At clos- 
ing velocities of 800 and 1,000 feet per second, one 
or more possible rendezvous trajectories achieved 
a positive value of range rate somewhere during 
the trajectory. This condition means that, for 
a while, the ferry vehicle was moving away from 
the space station. Such a phenomenon occurs 
when the ferry vehicle occupies a position ahead 
of and above the space station and rendezvous 
occurs on the descent from apogee. This condition 
leads to a largo negative value for a . (See cases 
illustrated in fig. 1(b).) If a display of range rate 
against range were the only information supplied 
to a human pilot during the midcourse phase of 
rendezvous, such a trajectory as this might lead 
to some consternation and doubt that rendezvous 
would indeed occur. 



(a) AF=400 feet per second. 

Figure 10. — Variation of range rate with range for a ferry 
vehicle during rendezvous with a station in a 300-mile 
circular orbit. All trajectories are in the orbital plane 
of the station 03=0). 



(b) AV=600 feet per second. 
Figure 10. — Continued. 


In all cases it may be noted that range rate 
was essentially constant during the last few miles 
before rendezvous. This condition, which is 
synonymous with a collision course, is found to be 
of interest in the terminal phase of a rendezvous. 
Such a condition is illustrated in figure 10(e) for 
the case where the closing velocity was 600 feet 
per second. Shown in figure 10(e) is an enlarge- 
ment of the last 60 miles for three of the trajec- 
tories shown in figure 10(b). It may be seen that 
during the last 60 miles the largest variation in 
range rate for any of the trajectories was 125 feet 
per second. In the last 12 miles (at a closing 
velocity of 600 feet per second, this value represents 
about 100 seconds of time) changes in range rate 
are virtually undetectable. Although the distance 
at which this condition occurred varied with the 
rate of closure, this trend is essentially the same 
whenever the ferry vehicle is on a collision course 
with the station. 
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(c) AF=800 feet per second. 
Figure 10. — Continued. 


Out-of-plane trajectories. — A number of the 
trajectories described in the previous section were 
also calculated with an initial out-of-plane velocity. 
For an in-plane velocity of 400 feet per second and 
a closing angle a of 0° trajectories were obtained 
with values of z at rendezvous of 200, 400, and 
600 feet per second. These values of i represent 
/3 angles of 26.6°, 45.0°, and 56.3°, respectively. 
These three cases are shown in figure 11(a) in 
plots of y against x and y against z. It may be 
seen that each 200 feet per second of lateral 
velocity represents a maximum change in out-of- 
plane distance of about 34 miles. It may also be 
seen that the variation of y with x is virtually 
independent of any variation in 2 . The approxi- 
mate equations of motion derived in appendix B 
resulted in solutions where 2 was independent of 
x and?/. The results of figure 11(a) indicate that 
the approximation (that z is independent of x and 
y) is fairly accurate, at least for half an orbital 
period. As a result of this condition, equation 

597870 — '61 3 
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(d) A V= 1,000 feet per second. 
Figure 10. — Continued. 



(e) Enlargement of three trajectories during the last 60 
miles prior to rendezvous. AU=600 feet per second; 
0 = 0 . 


Figure 10. — Concluded. 
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y 



(a) A V cos 0 — 400 feet per second ; a — 0. 

Figure 11. — Variation of out-of-plane trajectories in terms of x and y and z and y for the case of a ferry during rendezvous 

with a station in a 300-mile circular orbit. 


(B8) of appendix B may be used to compute a 
relationship for the maximum out -of -plane dis- 
tance obtained. For a given orbital angular veloc- 
ity of the station w and a given out-of-plane 
velocity z 0 evaluated at the time the ferry passes 
through the orbital plane, the maximum out-of- 
plane distance that can be obtained is 


2 max 


( 1 ) 


The physical analogy here is the simple harmonic 
motion of an undamped mass-spring system having 
a natural frequency of a>. The motion takes place 
normal to the plane of the orbit and the total 
energy of this motion is conserved. Thus, the 
maximum kinetic energy (z is a maximum when z 
is zero) is equal to the maximum potential energy 
( z is a maximum when i is zero). 

In the second example considered, the in-plane 
velocity was 600 feet per second and the out-of- 


plane velocity was 200 feet per second. These 
trajectories are shown in figure 11(b) for values of 
a of 30°, 0°, and —30°. Similar trajectories are 
shown in figure 11(c) where the out-of-plane 
velocity was 400 feet per second. A case was 
also calculated where the in-plane velocity was 
800 feet per second and the out-of-plane velocity 
was 200 feet per second and this case is shown in 
figure 11(d). It may be seen in each of these 
cases that the aforementioned approximate rela- 
tionship between z maz and z 0 is valid for the selected 
range of values of the in-plane component of the 
closing velocity AV cos 0. For each of these 
trajectories the angular distance traveled by the 
station, the launch velocity, and flight-path angle 
have been recorded in the figure. The reader 
should be cautioned, however, against a close com- 
parison of these numbers with those given in 
figure 9, since in both instances the numbers were 
obtained by interpolating between numerical data 
points. 



4* deg 

r u deg 

V L , ft/se< 

93.78 

2.33 

25,980.7 

99.84 

2.62 

25,903.2 

26.55 

2.20 

25,980.7 


100 



0 



0) 



1 -100 

- 





-200 

- 


1 -300 



300 -100 

^T r cos 

0=600 feet per second; AVsin 0 = 

E 11.— 

Continued. 


4> deg 

> 4 , deg 

V^ % ft/se< 

94.07 

2.19 

25,982.9 

00. 13 

2.60 

25,905.5 

27.00 

2.18 

25,983.0 


100 

n 

r~ 


u 

Ui 

0> 



i -loo 

- 







18 


TECHNICAL REPOET R-87— NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


a, deg 
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0 
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.10 26,221.5 



(d) A V= 824.6 feet per second; 0=14.0° (AV cos 0=800 feet per second; A Vain 0=200 feet per second). 

Figure 11. — Concluded. 


Boundaries of launch and rendezvous con- 
ditions. — In the discussion of the foregoing 
trajectories, it was noted that, under certain 
launch conditions and along certain trajectories, 
rendezvous could be achieved with a station in a 
300-mile circular orbit. With the orbit of the 
station and the launch altitude specified, the con- 
ditions at launch and at rendezvous will lie within 
certain definite boundaries. There are a number 
of ways of plotting such boundaries. If one con- 
siders the in-plane and out-of-plane motions as 
being independent and their effects additive, the 
boundaries of launch and rendezvous conditions 
may be plotted as a function of the angle of closure 
a. Also of interest is the range of launch velocity 
and flight-path angle commensurate with a given 
closing velocity. Plots illustrating these varia- 
tions are given in figure 12. 

If a limiting value of in-plane closing velocity 
AV cos ft is specified, the range of conditions at 
launch ( V L and y L ) which will achieve rendez- 
vous without exceeding the selected limiting value 
of closing velocity may be obtained from figure 
12(a). The region of possible launch conditions 
is represented by a half-crescent-shaped area de- 
fined by the three lines: the specified maximum 
allowable value of AV cos 0, the axis 7z,=0°, and 
the line a— 0°. Also plotted in figure 12(a) are 
lines showing the angle of closure a corresponding 
to the various combinations of permissible launch 
conditions. Positive values of a result when 


Al/ cos 0, ft/sec 



(a) Variation of the launch conditions showing lines of 
constant in-plane closing velocity and lines of constant 
in-plane closing angle. 

Figure 12. — Boundaries of ferry launch and rendezvous 
conditions for the case of a station in a 300-mile circular 
orbit. Launch refers to the condition of boost burnout, 
which is assumed to be 60 miles above the surface of the 
earth. 

rendezvous is achieved in the ascending phase of 
the transfer trajectory (see fig. 1); negative values 
are associated with rendezvous in the descending 
phase. 


A STUDY OF THE POSITIONS AND VELOCITIES OF A SPACE STATION AND A FERRY VEHICLE 


19 


The line a=0° is not a physical boundary but 
rather a line of symmetry. The two plrysical 
boundaries are lines of constant AV cos ft and the 
line of y L — 0°. The maximum allowable AY 
cos (3 is a measure of the maximum fuel carried by 
the ferry for injection into orbit. The region 
Yi<0° is excluded since a ferry launched at 60 
miles with a sufficiently small negative flight-path 
angle would reach a perigee within the atmosphere 
and again pass through the 60-mile altitude with 
a positive flight-path angle. It is this positive 
flight-path angle and the associated velocity which 
must fall within the specified boundary. (If the 
launch altitude were sufficiently large so that drag 
could be ignored, the magnitudes of y and of V 
would be unchanged and yz=0° could be treated 
as a line of symmetry. The launch times, how- 
ever, would be different.) 

Figure 12(a) also illustrates that, if one specifies 
the (circular) orbit of the station and then specifies 
two of the variables, for example, V L and y L , the 
trajectory conditions of the ferry as it passes 
through the orbital altitude of the station are fixed. 
(The quantities AV and a will be significant only 
if the station happens to be at the intersection of 
the orbits at that time.) Because of this restric- 
tion, the relative position of the two bodies at the 
time of launch is a dependent variable and cannot 
be independently chosen. Thus, there are only 
two independent variables to the problem once 
the station’s orbit is established. However, those 
two variables do not have to be V L and y L . For 
launch guidance one may prefer to specify a maxi- 
mum value of AV and a measured position of the 
ferry and then adjust the burnout conditions of 
the boost vehicle to satisfy the rendezvous 
conditions. 

The effects of time on the launch conditions are 
covered in a subsequent section. Although figure 
12(a) shows the necessary launch and correspond- 
ing rendezvous conditions jointly, additional in- 
sight into the boundaries may be obtained by cross 
plotting several of the variables. (See figs. 12(b) 
and 12(c).) 

In figure 12(b), the launch velocity and flight- 
path angle are plotted as a function of the angle of 
closure a. The boundaries of acceptable launch 
conditions lie between the line 7 z,= 0 ° and the 
maximum allowable value of AV cos /3. Only half 
the boundaries are shown since is plotteff on 
the abscissa. The contours indicate that, for a 




(b) Boundaries of inertial velocity and flight-path angle of 

ferry at time of launch as required for rendezvous. 

Figure 12. — Continued. 

given value of the in-plane closing velocity AV 
cos V L is a minimum and y L is a maximum at 
a=0°. As | a: | is increased, V L increases and ]y L 
decreases. 

In figure 12(c) the relative in-plane closing 
velocity and the ferry flight-path angle y f at ren- 
dezvous are plotted as functions of a. Again the 
appropriate boundaries or contours are marked. 
The variation of AV cos 0 with a is symmetric 
about the line a— 0° while the variation of y f 
with a is antisymmetric (positive y f with positive 
a). These boundaries show that, for a relative 
in-plane closing velocity of 400 feet per second, 
the angle of closure must lie between ±25° and 
the \y f \ must lie in a very narrow corridor between 
0.38° and 1.0°. At a relative closing velocity of 
1,000 feet per second, the relative angle of closure 
may be as much as +78° and the flight-path angle 
may be as much as ±2.23°. 
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(c) Boundaries of the relative closing velocity and flight- 
path angle of the ferry at the time of rendezvous. 

Figure 12. — Concluded. 

The boundaries given in figure 12 were ob- 
tained from two sources: the trajectories computed 
by the equations of motion of appendix A and the 
well-known closed-form solutions to Kepler's 
equations. The results obtained by the two 
methods were checked, one against the other, for 
possible errors. Since similar boundaries for any 
specified orbit may be computed with fairly simple 
equations, these equations and the considerations 
employed in their development are given in 
appendix D. It should be noted that these calcu- 
lated boundaries do not involve the out-of-plane 
component of the closing velocity nor do they 


involve the position and velocity as a function of 
time. For instance, the relative positions of the 
station and the ferry at any given time cannot be 
calculated in closed form. 

Variation of conditions with launch time. — 
For any given position of the space station with 
respect to the ferry at the time of booster burnout, 
the inertial velocity and flight-path angle of the 
ferry become uniquely specified if rendezvous or 
near rendezvous is to be accomplished with only 
small midcourse velocity corrections. For each 
second of delay in the time of firing of the booster 
on the ground, the station in a near-earth orbit 
will change its position by about 5 miles. There- 
fore, it appears certain that the conditions at the 
time of booster burnout must be varied to com- 
pensate for any delay in time. Such corrections 
could be fed into the booster guidance system 
right up to the instant of ground firing. Figure 13 
shows the manner in which these conditions might 
vary for the case where the station is in a 300-mile 
circular orbit. Shown in figure 13(a) is the 
variation of ferry launch velocity and flight-path 
angle with respect to the angular separation 
between the space station and the ferry for differ- 
ent values of the in-plane closing velocity, AT" 
cos 0. The separation angle A0 is positive when 
the station is ahead of the ferry. (See appendix 
C.) The abscissa is also marked in terms of r e A0 
which is the separation distance between the ferry 
and the space station measured in statute miles 
on the surface of the earth. It may be seen that, 
with each increase in AV cos j3 (over the range 
investigated), somewhat larger values of separa- 
tion distance or delayed times in launch may be 
accounted for. Also, if launch is attempted at 
the earliest possible time and the in-plane closing 
velocity is limited to 800 feet per second, the 
separation distance can decrease from 500 miles 
to 50 miles, and thus delays on the order of 90 
seconds may be accounted for. 

Depending upon when the vehicle is launched, 
rendezvous may occur from one-quarter to three- 
quarters around the earth from the launch posi- 
tion. The plots of figure 13(b) show the variation 
in the distance traveled around the earth from 
launch to rendezvous for certain specified rendez- 
vous conditions. In figure 13(b) are the values 
of velocity and flight-path angle of the ferry 
plotted against 0 L (that is, the position of the 
station at the time of booster burnout). Since at 



A STUDY OF THE POSITIONS AND VELOCITIES OF A SPACE STATION AND A FERRY VEHICLE 21 



r e AQ, miles 

(a) Relative position between ferry and station at launch. 

Figure 13. — Variation of the launch velocity and flight- 
path angle of the ferry as a function of the relative and 
absolute positions of the station and ferry at the time of 
booster burnout. The variations are shown for several 
values of the closing velocity. 

rendezvous 0=0°, negative values of 9 /y are 
shown. Contours have been drawn for launch 
conditions leading to rendezvous at specified 
values of the closing velocity AV cos (3. As 
previously shown in figure 12, for each closing 
velocity there is a minimum launch velocity which 
is associated with the angle of closure 0°. 
Contours through these minimum points have 
been drawn and are indicated on the figure as 
a— 0°. The solid line, marked a=0°, gives the 
value of 0 of the station at the time of ferry 



0^, deg 


(b) Absolute position of station at launch. 

Figure 13.— Concluded. 

launch. The figure shows, for example, that, 
with a closing velocity of 1,000 feet per second, 
the station is intercepted anywhere from 63° to 
246° around the earth from its position at ferry 
launch depending on the launch conditions. For 
a minimum launch velocity of 25,525 feet per 
second, rendezvous will occur about 73° around 
the earth from the launch point and at rendezvous 
a=0°. 

STATION IN AN ELLIPTIC ORBIT 

Trajectories.— A limited number of rendezvous 
trajectories were calculated for a station in an 
elliptic orbit. In order to have some basis of 
comparison with a 300-mile circular orbit, the 
elliptic orbit was chosen to have a perigee at 100 
miles and an apogee at 500 miles. The two orbits 
are illustrated in figure 14. 

Rendezvous was investigated at four positions 
along the elliptic orbit: at perigee (0—0°), at 
apogee (0=180°), and at the intersection of the 
orbit with the latus rectum (0=90°, 270°). 
These rendezvous positions are also illustrated 
in figure 14. It was assumed that at rendezvous 
the out-of-plane velocity was 400 feet per second 
and the in-plane velocity was 600 feet per second 
for a total AV of 721 feet per second. For all 
cases then, 0 was 33.69°. The in-plane angle of 
closure a was varied from zero in both negative 


y y miles 
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Rendezvous 


Elliptic--.^ 


0 = 180 ° 


Return 


0 n 90 ° 



0 = 270 ° 

Figure 14. — Circular and elliptic orbits. 


and positive directions until (in general) a value 
was reached where a launch from an altitude of 
60 miles was not possible. The results of these 
calculations are shown in figures 15 to 18. 

The inertial velocity of the station at the peri- 
gee of its orbit is 26,220.47 feet per second. The 
apogee velocity of a ferry on a trajectory whose 
perigee is 60 miles and whose apogee is 100 miles 
above the earth is 25,562.25 feet per second. 
Thus, the minimum closing velocity of the ferry 
and station at the station’s perigee is 658.22 feet 
per second. Therefore, when rendezvous at the 
station’s perigee with a closing velocity of 600 
feet per second was considered, none of the tra- 
jectories satisfied the required launch conditions 
A/=60 miles. However, a number of trajectories 
are shown in figure 15 for cases where a varied from 
50° to —50°, Shown on the (a), (b), and (c) parts 
of figure 15 are the y,x variation, the y,z variation, 
and the R,R variation, respectively. Listed in a 
table in the figure is the minimum altitude reached 


l , deg 
50 
40 
30 


68.69 
68.99 
71 .94 


0, deg 
316.09 
307.38 
290.27 


V fy ft /sec 

26 . 038.7 
25 , 962.2 

25 . 882.8 



x , miles 


z , miles 


(a) In-plane variation. (b) Out-of-plane variation. 

Figure 15. — Trajectories of a ferry vehicle*during rendezvous with a station in a 100- to 500- mile elliptic orbit. 
Rendezvous occurs at 0 = 0° (perigee). AV cos 0 = 600 feet per second. 




A STUDY OF THE POSITIONS AND VELOCITIES OF A SPACE STATION AND A FERRY VEHICLE 23 



R \ miles 

ate with range. 

Figure 15. — Concluded. 

by the ferry h ftmin (perigee), the 0 position of 
the station, and the inertial velocity of the ferry 
at this condition. Each of these trajectories was 
terminated when the ferry was at the minimum 
altitude. In order to facilitate a study of these 
trajectories, most trajectories have been marked 
by a symbol when the station is at 270°, 180°, 
and 90°. Whenever the trajectories could not be 
marked clearly, no mark was used. 

In figure 16 are shown the trajectories obtained 
for the ferry vehicle when rendezvous occurred at 
0=90°. All the trajectories shown are possible 
rendezvous trajectories and ticks have been put 
at the end of each trajectory to denote this fact. 
The conditions at launch for each trajectory are 
tabulated in the figure. It may be seen that a 
varied from —60° to 100°. The lack of S 3 unmetry 
in the range of a is due to the fact that the inertial 
velocity of the station is not alined with the AT-axis 
of the station at the time of rendezvous. 

Figure 17 shows all the trajectories obtained 
when rendezvous occurred at 0=180°. At 
apogee, the inertial velocity of the station was 

23.868.79 feet per second while a ferry on a 60- to 
500-mile orbit will have an inertial velocity of 

23.806.80 feet per second at this point. Thus, 
a minimum relative closing velocity of about 62 


feet per second is theoretically possible. For the 
specified in-plane closing velocity of 600 feet per 
second, the spread in the available a range for 
rendezvous is somewhere between ± 80 ° and ±90°. 
Symmetry occurs because the inertial velocity of 
the station and the AT-axis are again colinear at 
rendezvous. The launch conditions are tabulated 
in the figure for each trajectory shown. Symbols 
are placed on most of the trajectories to indicate 
the point at which 0 was 90°. 

In figure 18 the trajectories which led to ren- 
dezvous at 0=270° are shown. The trajectories 
are again marked in the same fashion as shown in 
figures 16 and 17. Symbols are used to indicate 
the position of the ferry when the station was at 
0 of 180° and 90°. In this case a varied from 
60° to -100°. 

It is interesting to note that the maximum 
available spread in a is almost identical for ren- 
dezvous at 0 of 90°, 180°, and 270°. In each 
of these cases a varied through about 160°. This 
result was surprising in that it was thought that 
the spread in a at 0=180° would be notably 
larger than that at 0 of 90° and 270°. 

Considered individually, the trajectories leading 
to rendezvous with a station in an elliptic orbit 
are very similar to those obtained when the sta- 
tion was in a circular orbit. However, if con- 
sidered as families of trajectories, the families are 
markedly different. One difference is primarily 
due to the fact that both the altitude and angular 
velocity of the rotating coordinate system are 
going through a nearly sinusoidal, rather than a 
constant , variation with time. In figures 16, 17, 
and 18, it may be noted that, if a line were drawn 
connecting the launch points of each of the tra- 
jectories, it would have a cycloidal character. 
With the station in a circular orbit this locus of 
launch positions was an arc of a circle of radius 
4,020 statute miles. (See fig. 9.) 

The 400-feet-per-second out-of-plane velocity 
component again produced a lateral displacement 
of the ferry of about 70 ±3 miles. It should be 
noted, however, that, when rendezvous occurred 
at 0=90°, the launch position was less than a 
quarter of the way around the earth and the lateral 
range of less than 70 miles is obtained. This con- 
dition implies that for these cases the out-of-plane 
velocity 7 of more than 400 feet per second must be 
used if a lateral displacement of 70 miles is 
required. 
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(a) In-plane variation. (c) Variation of range rate with range. 

(b) Out-of-plane variation. 

Figure 16. — Trajectories of a ferry vehicle during rendezvous with a station in a 100- to 500-mile elliptic orbit. Rendez 

vous occurs at 9 = 90°: AV cos 0 = 600 feet per second. 
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Boundaries of launch and rendezvous con- 
ditions, — The results of figures 15 to 18 and the 
equations of appendix D have been used to sum- 
marize boundaries of launch and rendezvous 
conditions for the case where a station is in a 
100- to 500-mile elliptic orbit. These boundaries 
are shown in figure 19. In figure 19(a) the range 
of possible launch conditions is shown for ren- 
dezvous with an in-planc velocity component of 
600 feet per second. Also shown in figure 19(a) 
are the variations of launch conditions when the 
closing velocity is varied but the angle of closure 
a is held constant at zero. Note that for rendez- 
vous at 0 = 90° and 270° that the a = 0° line lies 
within the boundary. 

In figure 19(b) the boundaries of the conditions 
at launch are shown in terms of V L and y L of the 
ferry plotted against the angle of closure a. In 
figure 19(c) the boundaries of rendezvous condi- 
tions, AV cos 0 and y / for the ferry are shown. In 
all but one case the boundaries are shown for 
rendezvous at 0—90°, 180° and 270° only. The 
one exception is in figure 19(c) in the variation of 
AV cos 0 with a for the case of 0=0° (or 360°). 
It may be seen that this contour does not intersect 
the line for AV cos 0=600 feet per second. 

Parallel to the presentation made in figure 10, 
the boundaries of figure 19 are given in terms of the 
in-plane closing velocity with the assumption that 



(a) Variations of the launch conditions showing lines of 
<*=0° and AV cos 0=600 feet per second. 

Figure 19. — Envelopes of ferry launch and rendezvous 
conditions for case of the station in a 100- to 500-mile 
elliptic orbit. Envelopes are shown for rendezvous at 
0 = 90°, 180°, and 270° for in-plane closing velocities 
of 600 feet per second or less and for ti. =0°- 



Boundary © at rendezvous, deq 

AV cos 0 = 600 ft/sec 90,180,270 



(b) Boundaries of launch conditions as a function of the 
angle of closure a. 

Figure 19. — Continued. 

the results were not dependent on t lie out-of-plane 
velocity. With the station in a circular orbit, 
this assumption was justified by showing that the 
in-plane and out-of-plane motions were virtually 
independent and could be treated separately. 
This condition has also been verified in the present 
case. A comparison was made of the results taken 
from figure 18 (with an out-of-plane velocity of 
400 feet per second) with those calculated with a 
zero out-of-plane velocity and with boundaries 
calculated by the method of appendix D. (In- 
plane considerations only were used.) Virtually 
no difference could be detected in the launch condi- 
tions. This fact can be appreciated by noting that 
the vector addition of a 400-foot-per-second ve- 
locity component normal to a 25,000-foot-per- 
second component produces an increase in velocity 
of about 3 feet per second. Thus, within the 
reading accuracy the plots of figure 19 may be 
applied to out-of-plane trajectories, but caution 
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(c) Boundaries of rendezvous conditions as a function of 
the angle of closure a. 

Figure 19. — Concluded. 


should be used if /3)>6 0°. 

In figure 19, the acceptable range of variables 
again lies between boundaries defined by the con- 
tours for y L — 0° and the maximum acceptable 
(in-plane) closing velocity, arbitrarily taken to be 
600 feet per second. In figure 19(c) the 7 ^= 0 ° 
contours have been extended out to the condition 
where AV cos ft— 1,000 feet per second in order to 
show the extreme sensitivity of this contour at the 
larger values of closing velocity. In figure 19(b) 
the hatching has been left off the AV cos ft bound- 
ary of launch velocities for the sake of clarity. 

Examination of the variations in these bounda- 
ries for rendezvous at the various positions of the 
station indicates that in the sector of the ellipse 
defined by 270° ^ 0 ^ 90° (the apogee sector) : 

(1) The minimum values of launch and closing 
velocity do not vary greatly as 0 changes. These 
minimum values occur between values of a of 
± 20 °. 

(2) Kelative velocity and flight-path angle at 
rendezvous change very little with a over a range 
of about 100° (fig. 19(c)). At the extreme end of 


the a range, however, the variations of AV cos ft 
and 7/ change very rapidly with small changes in a. 

(3) For a AV cos 0^200 feet per second, the 
total range of a available for rendezvous for any 
given value of closing velocity is about the same 
for any value of 0. 

(4) The boundary of launch conditions (fig. 
19(a)) is very narrow for rendezvous at 0=180° 
but opens up considerably for rendezvous at 
0=90° and 270°. 

(5) With a maximum in-plane closing velocity 
component of 600 feet per second, the launch 
velocity may be reduced by as much as 500 feet 
per second below that required for a minimum 
closing velocity rendezvous. This result is similar 
to that noted for the case of the station in a circular 
orbit. 

Very little advantage can be seen for rendezvous 
in the sector between 27O°<0<9O° (the perigee 
sector) since conditions are changing rapidly 
along that part of the orbit. It can, for example, 
be seen in figure 19(c) that the minimum in- 
plane closing velocity changes from 108 feet per 
second to 660 feet per second between the in- 
tersection with the latus rectum and perigee. 
However, it is interesting to note that the 7 ^= 0 ° 
boundary for rendezvous at perigee also exhibits 
the weak variation of AV cos ft with a displayed 
by the curves for rendezvous at the other three 
positions. 

Variation of conditions with launch time. — 

In the case of the station in a circular orbit, the 
variation of the launch conditions and the time 
available for launch were limited only by the 
desired closing velocity. The down-range posi- 
tion at which rendezvous was to occur was not 
important because there was only one possible 
orientation of the orbit of the station with re- 
spect to the circular earth. With the station 
in an elliptic orbit this limitation no longer exists. 
Because the station velocity and altitude are 
changing continuously, the orientation of the 
station’s orbit and down-range position of the 
station at which rendezvous is to occur will 
directly affect the launch conditions. Although 
each position of the station is unique, the trend 
of the change in launch conditions with the value 
of 0 at rendezvous can be seen in figure 20. 

Figure 20(a) shows the relative position of the 
station and ferry at the time of ferry launch. 
The scales shown along the abscissa give both 


i 


, ft/sec 
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(a) Boundaries of launch conditions as a function of the 
relative position of the ferry and station at launch, A0 
and tvAB. 

Figure 20.— Variation of launch velocity and flight-path 
angle for the ferry as a function of the absolute and 
relative positions of the station at launch. Contours 
represent a maximum in-plane closing velocity of 600 feet 
per second. 

the separation angle measured from the center 
of the earth A0 and the separation distance 
measured oyer the surface of the earth r e A9. 
Positive values of A0 are used when the station 
is ahead of the ferry at the time of launch. Fig- 
ure 20(b) shows the absolute position of the station 
at the time of ferry launch when rendezvous is 
to occur at 600 feet per second at 0 = 90°, 180°, 
and 270°. The values of 0^ shown on the abscissa 
are in agreement with the definition of 0 for 
an elliptic orbit; namely, 0z,=O° means that the 
ferry launch occurred when the station was at 
the perigee of its orbit. It may be noted that 



(b) Boundaries of launch conditions as a function of the 
absolute position of the station at launch 0 l. 

Figure 20. — Concluded. 

the curves of figure 20 arc the contours for AV = 
600 feet per second of figure 19 and that the 
y L ~0° contours could be used to close the 
boundaries. 

Figure 20(a) indicates that launch times are 
very limited if rendezvous is to be made at 0—90°. 
For this case the station must be somewhere 
between 125 statute miles behind to 130 statute 
miles ahead of the ferry at the time of launch. 
This total spread of 255 statute miles is compa- 
rable to roughly 380 miles when the station is 
in a circular orbit, where in both cases the in- 
plane closing velocity is 600 feet per second. 
Rendezvous at 0 = 180° and 270° offer a larger 
spread in launch time. 

In figure 13(a) it was shown that, with the 
station in a 300-mile circular orbit, most of the 
possible launch conditions occurred with the 
ferry behind the station. For the station in 
a 100- to 500-mile elliptic orbit, figure 20(a) 
shows that the relative positions at launch are 
about equally divided between the ferry ahead 
of or behind the station. 

Figure 20(b) shows that a rendezvous at 0=90° 
can be obtained by launching from a position 
between 0= —41° and 44° and that the launch ve- 
locity changes less than 300 feet per second between 
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the two extremes. Thus a rendezvous could 
occur at a position 46° (90° — 44°) down range 
from the launch if all other conditions were met. 
Rendezvous at 180° and 270° positions require 
progressively longer down-range distances. 

DISCUSSION 

In the foregoing sections, typical trajectories 
leading to rendezvous with a station in an elliptic 
and a circular orbit were illustrated. Boundaries 
which define the allowable launch and rendezvous 
conditions were also developed. There is, of 
course, nothing unique about the particular 
orbits chosen as examples. The range of allowable 
near-earth orbits is usually restricted to the 
region above the sensible atmosphere of the earth 
and below the Van Allen radiation belt. An 
altitude of 300 miles has often been considered 
a mean value for this region. 

The trajectories were limited to those which 
produced a closing velocity at rendezvous of 1,000 
feet per second or less. This limitation was im- 
posed since, according to the formula for an 
impulsive velocity change 

^=l-exp(-^£)«^ (2) 

a AV of 1,000 feet per second represents about 
12.5 percent of the payload of the ferry ( I sp —25Q 
seconds). The results show, however, that this 
increase in the closing velocity is not entirely lost, 
since in most cases there is a reduction in the 
launch velocity. If the closing velocity is in- 
creased from the minimum of 360 feet per second 
to a value of 1,000 feet per second, the launch 
velocity may be reduced from a value of 26,125 
feet per second to a value of 25,525 feet per second. 
(See fig. 10(b).) Thus, an increase in the relative 
velocity at rendezvous of 640 feet per second may 
be compensated by as much as a 600-foot-per- 
second reduction in launch velocity. This large 
reduction in launch velocity is, of course, realized 
only at one instant of launch time and any error 
in this launch time will lead to a requirement for 
somewhat larger launch velocities for the same 
closing velocity. 

With reference to launch times, it is interesting 
to note in figure 13(a) that, regardless of the clos- 
ing velocity, the station cannot be more than 
about 540 miles ahead of the ferry at the time of 
launch in order to achieve rendezvous in a 300- 


mile circular orbit. Increasing the closing velocity 
allows the ferry to be launched at an earlier time 
but does not allow the ferry to be launched at a 
later time. 

The results show that the out-of-plane motion 
of a ferry may be treated as virtually independent 
of the in-plane motions. If rendezvous is to 
occur after the ferry has passed more than a 
quarter of the way around the earth, the ferry 
may be launched in a direction parallel to the 
orbital plane of the station. At a point a quarter 
of the way around the earth, the ferry will pass 
through the plane of the station’s orbit and an 
attempt should be made to bring the out-of-plane 
velocity i to zero. The cost of this maneuver will 
be a function of the out-of-plane distance at 
launch and is given by equation (1). If rendezvous 
is to occur at a point less than a quarter of the way 
around the earth from launch, the ferry must have 
a z component at launch sufficient to insure that 
the ferry will intersect the plane of the station’s 
orbit before or at the time of rendezvous. One 
reason for this prior intersection, other than sim- 
plicity of injection guidance, is found in the fact 
that the inertial velocity of the ferry is lower 
before the ferry is injected into orbit than it is 
afterwards. Therefore, the out-of-plane velocity 
corrections will be smaller. 

In the case where the station is in an elliptic 
orbit, it appears that the variation in the launch 
conditions could become very critical. The results 
of the one elliptic orbit considered indicate that 
perigee would be the least desirable position to 
attempt rendezvous and near apogee, the most 
desirable. Near apogee, the available launch time 
is relatively large, the required injection velocity 
is a minimum, and the conditions at rendezvous 
are relatively invariant with the angle of closure 
over a fairly wide range of conditions. Rendezvous 
between perigee (0=0°) and the ascending inter- 
section with the latus rectum (0 = 90°) offers one 
advantage in the form of having the rendezvous 
occur relatively close to the launch site (on the 
order of 45° down range compared with 90° to 
180° around the earth for rendezvous at apogee). 

The choice of positions of the station at the time 
of rendezvous may, of course, be somewhat aca- 
demic since in practice the orbit of the station 
must also pass relatively close (laterally) to the 
launch site at the time of launch. Thus, wdien 
the station is in an elliptic orbit, some compromise 
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may be necessary between an orbital pass that lias 
favorable lateral characteristics and one that has 
favorable orientation characteristics. With a cir- 
cular orbit, this orientation problem does not arise. 

In programing the launch conditions to com- 
pensate for delays in launch time, rendezvous with 
a station in an elliptic orbit appears to be more 
difficult than that of a station in a circular orbit 
since the altitude, velocity, and down-range posi- 
tion of the station at which rendezvous is to occur 
arc constantly changing with time. The results 
shown in figure 20 indicate that these changes arc 
not severe for the example orbit but would prob- 
ably become more severe as the eccentricity of the 
station's orbit is increased. 

CONCLUSIONS 

It may be concluded that, for any orbit that a 
space station may have, there is a family of trajec- 
tories leaving the earth's atmosphere with different 
velocities and along different flight paths which, 
without further thrusting, will terminate at some 
instantaneous position of the station. These fam- 
ilies may be computed and the launch conditions 
required for rendezvous determined. The equa- 
tions and procedures have been developed for 
computing these trajectories and the proper 
launch conditions (at time of booster burnout). 

Typical calculations have been performed for a 
station in a circular orbit 300 miles above the earth 
and for a station in an elliptic orbit having a perigee 
100 miles and an apogee 500 miles above the earth. 
These numerical studies indicate that the launch 
velocity may be reduced as the closing velocity is 
increased but, in no case, is the sum of the two less 


than that obtained for the Hohmann ellipse. As 
the allowable closing velocity is increased, the 
range of possible launch and rendezvous conditions 
increases and the time available for launch in- 
creases. If the closing velocity between the two 
vehicles at rendezvous is limited to about 500 feet 
per second above the minimum, rendezvous will 
generally occur with the station overtaking the 
ferry. However, with the station in an elliptic 
orbit, rendezvous can occur with the ferry over- 
taking the station from above or below. In 
either case the rate of closure approaches a con- 
stant and remains essentially constant during the 
last few minutes before rendezvous. This condi- 
tion, which is indicative of a collision course, 
could provide a criterion for guidance during the 
terminal phase. 

If the station is in an elliptic orbit, rendezvous 
appears to be more desirable near apogee than 
near perigee. Furthermore, it appears that launch 
guidance for rendezvous with a vehicle in an 
elliptic orbit will be more critical that that for one 
in a circular orbit ; however, once launched, the 
range and variation of available rendezvous condi- 
tions make rendezvous with a station in an elliptic 
orbit more favorable if rendezvous is to occur 
somewhere near apogee. In either case, however, 
it appears that the prediction of desirable rendez- 
vous trajectories is feasible and that launch and 
probably midcourse guidance will bo necessary to 
allow for delays in launch time and to insure that 
the proper trajectory is established. 

Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., October 27, 1960. 
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APPENDIX A 


EQUATIONS OF MOTION OF A MASS AS MEASURED IN A ROTATING FRAME OF AXES WHICH FOLLOWS 

A KEPLERIAN TRAJECTORY 


The method of Lagrange is employed in deriving 
the equations of motion based on the kinetic and 
potential energy of the mass at any instant. 
Sketch 1 shows the coordinates and the velocity 





Sketch 1. 

notation employed. The position of the mass is 
given by 

r/=r s +R (Al) 


The velocity of the mass with respect to the 
center of the earth's gravitational field (taken as 
the absolute system) and with respect to the 
rotating x,y,z coordinate system is given by 

v -^-7ff+Tr+ axK < A2) 

where the prime denotes the change as seen from 
tiie origin of the coordinate system. The vectors 
in equation (A2) have the components 

D=(o,o,e) 

R ={x,y,s) 

r,= (x,y+r s ,2) 

~={-r s Q,r s , 0) 

where r s and 0 are not coordinates of the mass but 
are specified time-dependent parameters. The 
absolute velocity of the mass in terms of the coor- 
dinates x,y,z , and the parameters r SJ 0 is then 

\^i[x—(y+r s )Q]+\(y+r s +xQ)+k(z) (A3) 
The kinetic energy of the mass is then 

T=\m(y • V) 

x 2 +y 2 -\- z 2 +{(y+r s ) 2 +x 2 ]Q 2 +2yf s 
+r s 2 +2[x(r s +y)—(y+r s )x]ey (A4) 

The potential energy of the mass is given by 

U= — mgr f 

mg e r e * 

r; 

=—mg e r 2 [x 2 +(y+r s ) 2 +z 2 ]~ m (A5) 
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With the kinetic and potential energy thus 
specified, the Lagrangian 

L=T—U (A0) 

can be formed and the Lagrangian operation 


obtained : 


the semimajor axis: 



(A7) 

performed for the equations of motion in terms of 
the coordinates x, y , and 2 (which are q u q 2 , and 
q s in eq. (A7)). 

After the indicated operation of equation (A7) 
is made, the equations of motion are 


the eccentricity: 


^s.a + Ts, p 

and the angular momentum: 


C 9 =lgs*A( W )] 1/2 


(All) 


(A12) 


(A 13) 


X— (y+r s )Q-2(y + r,)Q-xQ 2 +g e r e 2 

T f III 

y+xQ+2xe+r s -(y+r s )&+g e r e 2 -^3 —=^ ► 
... , z T z 

S +9° r ‘7}=m J 

(AS) 

The symbols T x , T yy and T z denote the thrust 
forces acting on the mass in the £-, y-, and z- direc- 
tion of the station coordinate system. The 
radial distance from the center of the gravitation 
field to the mass is 

r f =[x 2 +(y+r s ) 2 +z 2 ] 112 (A9) 

The equations which specify the variation of r s 
and 0 are the Keplcrian equations: 


Equations (AlO) can then be replaced with 

t== k £ r ^ e)de 


r (e) = jd£kzf?L 

^ 1+e cos 0 


and the time derivatives 

• sin 0 

Vs ~~ A( 1-e 2 ) 

C s 2 6 COS 0 


0=^ 


A( \-e 2 ) r s 2 
C s 


0 = ■ 


2 C s 2 e sin 0 
m A( 1-c 2 ) r s 3 J 


(A 14) 


(A 15) 


r s 2 Q=C s 

where C s is the constant angular momentum of the 
space station. Since equations (AlO) can be 
solved in closed form, it is normally more con- 
venient to specify the conic section in terms of its 
characteristic constants and the initial position of 
the station. For an ellipse, one may define 
0—0° at the perigee and define the perigee r SfP , 
the apogee r s>a , and the initial position of the 
space station in the ellipse 0 O . With these 

definitions, the following relationships can be 


For calculating the line-of -sight range between the 
two vehicles and the time rate of change of that 
range, the expressions 

R = '\ / x 2 +y 2J rz 2 (Al6) 

n = B±>m±ii (A17) 

may be used. 

It should be noted that the velocities measured 
from the rotating frame of axes in the space station 
are not the same as those that would be measured 
from the absolute frame of axes fixed at the 
center of the earth. As can be seen from equation 
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(A3), the proper relationships are 

V x =x—(y+r s ) 0 

V„=y+r s +xQ 

V,= i 

v^vj+vj+v/. 


(Al8) 


where \=^ is given by equation (A3) and the 

vertical bars denote the absolute value of the 
resulting vector. The condition that the total 
energy remain constant requires that 


T-\~ U=*Ef= Constant 


(A20) 


where V Xi V Vj and V z are the components of the 
inertial velocity of the ferry in the direction of 
the X-, Y-, and Z- axes, respectively, as measured 
on the nonrolating earth. Velocity components as 
measured by a tracking station on a rotating earth 
could be obtained but these are not necessary to 
the subject of this paper. 

Two additional equations can be obtained which 
may prove useful if the differential equations are 
solved by numerical integration processes such as 
used by digital computers. Regardless of the 
orbit that is followed, if the forces are conservative 
(that is, aerodynamic forces and thrust are zero), 
the angular momentum and the total energy of 
the ferry vehicle must remain constant. The re- 
quirement that angular momentum be conserved 
requires that 

l r /X V|=6/=Constant (A19) 


Equations (A 19) and (A20) are the first integrals 
of the motion of the ferry vehicle in terms of the 
coordinates of the space station. If C f and E f are 
evaluated by using the initial conditions, at any 
point on the trajectory the degree of accuracy of 
the numerical solutions may be determined by 
again computing either C f or E f and comparing 
the percent of change. In the case of the angular 
momentum, it is suggested that the vector quanti- 
ties be evaluated in numerical form since the 
analytical expression is very cumbersome. How- 
ever, for the simplified case of coplanar motion of 
the ferry vehicle z—z— 0, the expression becomes 
simply 

Cf= [y + r s ] x— [y + r s ] x + [ x 2 + (y +r s ) 2 ]( 0 

(A21) 
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APPENDIX B 


EQUATIONS OF MOTION AND APPROXIMATE SOLUTIONS TO THE MOTION OF A MASS AS MEASURED 
IN A ROTATING FRAME OF AXES WHICH MOVES IN A CIRCULAR ORBIT 


For the case where the space station, and hence, 
the rotating frame of axes moves in a circular (or 
very nearly circular) orbit, the equations of ap- 
pendix A may be simplified. If 

r,— Constant 
0=Constant=co 



the differential equations (A8) become 

y+2ux+(y+r a )(^-a,^=^ (B2) 

2 + 2^f=-‘ 
rf m 

The radial distance to the ferry vehicle is defined as 

r f = [: x 2 + (y +r s ) 2 + z 2 ] 1/2 (B 3) 

and the orbital equations (AlO) give the expression 
for co in terms of the radius r s as 




(B4) 


In these equations it is an implicit assumption 
that the mass of the station, the mass of the ferry 
vehicle, the gravitational attraction between the 
two masses and any other planetary bodies are all 
negligible compared with the mass and gravitation 


field of the earth. If it is desired to go a step 
further and to consider cases where the ferry is 
never more than 100 or 200 miles from the station 
(in a circular orbit), very good approximate 
solutions can be obtained to equations (B2). The 
method given in reference 6 and first noted in 
reference 7 requires the expansion of rf z into a 
power series 


-1 = J_ 

*3 r 3 

f f • s 


[l -3+y+z'— 4V+,')- 



<B5) 


In the series, the terms of second order and higher 
are dropped. Substituting 

g e r 2 g e r 2 / ,, y 

r/ r s 3 \ r, 

« u 2 (l-3^) (B6) 


into equations (B2) and dropping the nonlinear 
terms containing ylr s gives a set of linear ordinary 
differential equations: 


x— 2c cy- 


m 


y+2ux-3a> 2 y=^ ► 

T 

Z + C0 2 Z = h 

m J 


(B7) 


On a nonthrusting trajectory, the position and velocity of the mass at any time are given in terms 
of the initial conditions by the homogeneous solutions: 

x=2 ( 2 --3 A sin cot- 2 ^ cos «<+( &y 0 -3 — ) coi+x 0 +2 & 

y=(% 32/„) cos w«+^ sin wt+ly 0 — 2 ^ (B8) 


z=z 0 cos cot+~ sin cot 
co 
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-=2 (2 — — tyo) cos wt+2 ^ sin at+6y„-3 - 

CO \ CO / CO *10 

-=—(2 — — 3u,^ sin cot-\-^- cos co£ ► 

0) \ o> / CO 


-=— cos cot—z 0 sin cot 
co co 


(B9) 


It may be noted that all these quantities, with the 
exception of x which increases linearly with time, 
are periodic in time. 

As a check on the numerical examples of this 
paper, the position and velocity of a mass (the 
ferry) ejected from the space station at time zero 
with the velocity components x 0 , y 0 , and z 0 can be 
obtained at any subsequent (or prior) times by 


setting x 0 =;?/ 0 = z o =0 in equations (B8) and (B9) 
and solving for x, y, z, x, y, and z at time t (or — <). 

On the other hand, if the position of the ferry is 
known at any time, for example, *= 0, the velocity 
components of the ferry relative to the space 
station necessary to rendezvous at some time r in 
the future can be obtained by putting equations 
(B8) in a slightly different form. Requiring 


that x=y~z ~ 0 when t~r yields 

’ 0 *1 pt sin cot — 3cor 2 — 2 cos cot 

2 cos cor— 2 sin cot 



] 


CO 

Vo 


K 


1 6cor — 6 sin cor” 


0 4—3 cos * 



z 

0—z o cos cor-1 — - sin cor 

CO 

Solving for the velocities in terms of the displace- 
ments gives 

x 0 _ x 0 sin cor+y 0 [6cor sin cor— 14 (1 — cos cor)] ^ 
co A 

Vo 2x 0 (1 — cose or)-\-y 0 (4 sin cor — 3cor cos cor) ^ 

co A 

%0 %0 

co tan cor ^ 


where 


A=3cor sin cor — 8(1 — cos cor) 


Again, it should be noted that the velocity com- 
ponents of the ferry as measured from the space 
station are not the same as the components meas- 
ured from a nonrotating earth. The proper 
relationships are obtained from equation (A3) as 


V x =x— (y+T s )co " 

V y =y+xco 

V'=i 

Vt^Vf+Vf+V*, 


(Bll) 


(BIO) 


APPENDIX C 

POSITION AND VELOCITY OF FERRY WITH RESPECT TO EARTH-FIXED COORDINATE SYSTEM 


The position of the ferry vehicle as seen from 
the center of the earth will be defined with the 
spherical coordinates r f , A0, and f as shown in 
sketch 2. The angle AO lies in the orbital plane 
of the station and the angle f is measured normal 


Y 



to the orbital plane. In terms of coordinates of 
the space station, these earth-measured coordi- 
nates are defined by 


AG = tan- 1 -~ 

y+r, 

r=tan ~ l [^+(2/+r 5 ) 2 ]^ 
r f =[x i + ( y +r s ) 2 + z 2 ] m 


(0 ^ A0<2jt) "1 


(— 7r<fgir) f 


(Cl) 


The velocity components of the ferry as meas- 
ured in a Cartesian frame of axes fixed at the 
center of the earth were given by equations (Al8) 
and (Bit) for the two classes of orbits considered. 
If these velocity components are rotated first 
through the angle AO and then through the 
angle f, the spherical velocity components are 
given by 

V T ~r f — (V x sin AO+Fy cos A0)"cos f+F 2 sin f"' 

F^ r/ ^-(F,sinAe+F,cosAO)sinf+F,cos f ► 

Ve~r f cos £AQ=V X cos AO— V y sin AO ^ 

(C2) 

and are shown in sketch 2. 

The velocity components Ff and F e define 
the plane of the local horizontal and the total 
tangential (horizontal) component of the velocity 

is 

(C3) 


and makes an angle 

77/^tan -1 


Fe 


(C4) 


with the plane of the station’s orbit. The local flight-path angle is defined as the angle between 
the resultant velocity vector and the local horizon and is given by 

7 /=tan _1 y- (Co) 

It should be noted that the in-plane angular position of the ferry with respect to the reference 
position 0=0 is given by the angle 0— A0=0 r . 
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APPENDIX D 

ANALYTICAL CALCULATION OF BOUNDARIES OF LAUNCH AND RENDEZVOUS CONDITIONS 


The boundaries of launch and rendezvous 
conditions, given in figure 12 for a circular orbit 
and in figure 19 for an elliptic orbit, can be calcu- 
lated analytically from relatively simple con- 
siderations. Since the procedure for calculating 
the boundaries is the same for either type of orbit, 
the equations are given for the more general case 
of the station in an elliptic orbit. The simplifica- 
tions obtained with the station in a circular orbit 
will then be noted. 

For a specified orbit of the station where 
rendezvous is to occur at the angular position 0, 
the conditions at that position are given by 


rendezvous are obtained at the condition r f =r s ; 
therefore, 


=2 g e r e 2 ( 


1 \ 

(D6) 

Zs 

r r.v+ r t,a) 

COS Y / = 

Zap 

r /.a 

Vf,v 

Vs. a 

(D7) 


The relative velocity components^are computed 
from the relations 

x=V s cos y s — V f cos y f 

y—V f sin y f —V s sin y s 


r s = 


A(l~e 2 ) 
l+€ COS 0 


Grid 


tan Ys— - 


€ sin 0 


9 1 + e cos 0 


(Dl) 

(D2) 

(D3) 


The expressions for A and e are given in appendix 
A. (See eqs. (All) and (A12).) 

In order to obtain the 7 ^ = 0° boundaries for 
rendezvous at the specified position of the station, 
the perigee of the ferry trajectory is held fixed at 
the launch altitude so that 

r /,p =r £=60+3, 960=4,020 miles 


The situation is shown in sketch 3. Note that 



only coplanar motion is considered here; there- 
fore, with 0=0, AV is used instead of AV cos 0. 
The velocity and angle of closure are then ob- 
tained from the relations 


and the apogee radius r f>a is varied from the radius 
of the station r $ over a range of values greater than 
r s . For each value of r SA that is chosen, the 
associated launch conditions are 


Vj} = V f ;?=‘2g e r e 2 (J— ^ ) 

V ' p i *f, a./ 


7^=0 


(D4) 

(D5) 


The velocity and flight-path angle of the ferry at 


A V 2 =x 2 Ay 2 

=TV+F/-2F,F /CO s (t/~Ys) (D8) 

tan a~% 
x 

F s cos Ys Fy cos 7/ 

Thus, equations (Dl) to (D3) give the station 
conditions at rendezvous; equations (D4) and (D5) 
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give the ferry launch conditions; equations (D6) 
and (D7) give the ferry rendezvous conditions; 
and equations (D8) and (D9) give the relative 
conditions at rendezvous. 

The constant AF boundaries are obtained by 
following the procedure in reverse. By holding 
AV constant and varying a, one obtains 

F/= V s 1 2 3 4 + AV 2 —2V s AV cos (a+y s ) (DIO) 


_V S sin 7 s +AF sin a 
tan y/ ~v s cos 7 S — AF cos a 


(Dll) 


Thus, AV and a are specified; equations (DIO) 
and (Dll) give the ferry conditions at rendezvous; 
and equations (Dl2) and (D13) give the ferry con- 
ditions at launch. This procedure may also be 
followed for contours of constant a and variable AV 

Extreme caution should be exercised in the use 
of equations (DlO) to (D13). When AF is close 
to its minimum possible value and also when y L 
approaches 0°, the equations become particularly 
sensitive to small computational errors. 

When the station’s orbit is circular, all the above 
equations still hold with the simplifications that 


Without calculating the perigee or apogee of the 
transfer ellipse, the conditions at ferry launch are 
obtained from the expressions of constant total 
energy and constant angular momentum: 


Ys~ 0 (D14) 

6=0 
t s ~A 


V^mV^-2 gjr* (L-£) (D12) 



COS 7 L = 


r s Vs cos 7 S — AF cos a 


(D13) 


and the absolute value of 0 at the time of ren- 
dezvous is no longer relevant. 
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